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ABSTRACT 



Aims. In this paper, we study a simple model of an orbiting protoplanet in a central magnetospheric cavity, the entry into such 
a cavity having been proposed as a mechanism for halting inward orbital migration. 

Methods. We have calculated the gravitational interaction of the protoplanet with the magnetosphere using a local model and 
determined the rate of evolution of the orbit. 

Results. The interaction is found to be determined by the outward flux of MHD waves and thus the possibility of the existence 
of such waves in the cavity is significant. 

Conclusions. The estimated orbital evolution rates due to gravitational and other interactions with the magnetosphere are 
unlikely to be significant during protoplanetary disk lifetimes. 
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1. Introduction 

The discovery of extrasolar giant planets orbiting close to 
their host stars ( Mayor & Queloz 1995; Marcy & Butler 
1995, 1998) with periods of a few days has led to the pro- 
posal that orbital migration during and post formation is 
responsible for them attaining short period orbits. This 
is because of the difficulties associated with forming such 
planets in situ, either in the critical core mass accumu- 
lation followed by gas accretion scenario, or the gravita- 
tional instability scenario for giant planet formation (see 
Papaloizou & Terquem 2006 and references therein for 
more discussion). However, estimates of disk driven mi- 
gration time scales have in general found these to be short 
compared to protostellar disk lifetimes (e.g. Nelson et al. 
2000). This in turn suggests the need for a mechanism to 
halt migration so preventing the protoplanet from falling 
into the star. Lin Bodenheimer & Richardson (1996) sug- 
gested that entry into a stellar magnetosphere close to the 
star could provide such a mechanism. This is through de- 
tachment from the disk. Indeed once the protoplanet is 
interior to the 2:1 resonance with the inner disk edge, in- 
teraction with the disk would be expected to cease through 
the lack of effective outer Lindlbad resonances (see e.g. Lin 
& Papaloizou 1993). 

However, there has as yet been no detailed analysis of 
the expected orbital evolution of protoplanets in magne- 
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tospheric cavities or demonstration of the effectiveness of 
this mechanism for halting inward orbital migration. This 
would require consideration of models for this inner region 
and a calculation of the interaction with the protoplanet 
and the effects on the protoplanet orbit. It is the purpose 
of this paper to make some first steps in this direction. 

Accretion on to T Tauri stars is indeed thought to be 
magnetically dominated close to the central star in some 
cases (Cameron & Campbell 1993; Menard et al. 2003). 
Surface field strengths of a few kilogauss (e.g. Safier 1998; 
Johns-Krull et al. 1999) have been estimated. According to 
models of accretion flows with their central regions domi- 
nated by the magnetic field of the central star ( e.g. Ghosh 
& Lamb 1978; Cameron & Campbell 1993; Matt & Pudritz 
2004) , at large radii the accretion flow is that of a normal 
viscous accretion disk ( e.g. Pringle 1981). The magnetic 
field penetrates the disk due to diffusion arising from the 
growth of various instabilities. At large radii the magnetic 
field is weak such that viscous stresses are more important 
than magnetic stresses. Then the effect of the magnetic 
field on the disk structure can be ignored. However, at 
smaller radii the dominance of magnetic stresses leads to 
the truncation of the disk at a radius, i?d, and the chan- 
neling of the flow to the star along stellar magnetic field 
lines. 

In order for accretion on to the central star to proceed 
the magnetic field must not disrupt the disk exterior to the 
corotation radius, where the angular velocity in the disk 
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coincides with that of the star, since the magnetic torque 
on the disk would impart angular momentum to the disk 
gas, making it unable to connect to the stellar field. If 
disk disruption occurs too far inside corotation the radius, 
there is an accretion torque acting to spin up the star. In 
general one can expect there to be an equilibrium state 
where spin up torques due to accretion and spin down 
torques due to magnetic stresses acting on the disk cancel. 
Because of the expected rapid increase of the magnetic 
stresses as the distance to the star decreases, one expects 
that i?d is slightly less than the corotation radius (Ghosh 
& Lamb 1978) when the system accretes onto the central 
star and is near equilibrium. 

This general picture of magnetic accretion has been 
found in recent three dimensional numerical simulations 
(see Romanova et al. 2004; 2006; Bouvier et al. 2006). 
It has been found that accretion along field lines takes 
place in funnel stream flows. When the stellar field is a 
dipole that is either not nearly aligned or strongly mis- 
aligned with the angular momentum axis, the magneto- 
spheric density is relatively low in the equatorial plane. 
For both small and large misalignment angles, relatively 
high density streams may occur there. The latter situation 
may be less favourable for stalling migration as interaction 
with the protoplanet is potentially stronger. 

The purpose of this paper is to consider a simple model 
of the stellar magnetosphere with a steady state accretion 
flow and to use it to calculate the protoplanet magneto- 
sphere interaction and the consequent protoplanet orbital 
evolution rate. This is to establish to what extent entry 
into such a magnetosphere is by itself sufficient to halt 
migration without the need to postulate additional effects 
such as residual coorbital torques resulting from the inner 
disk edge ( e.g. Masset et al. 2006). 

The plan of the paper is as follows: In section [2] we 
give the basic equations of the problem. In section [3] we 
consider perturbation of a steady state magnetosphere by 
a protoplanet in circular orbit with the aim of evaluating 
the exchange of energy with the orbit resulting from the 
perturbed gravitational forces. 

We carry out a purely local calculation of the proto- 
planet magnetosphere interaction. This confirms the im- 
portance of the existence of propagating waves in order to 
obtain a non zero energy exchange rate. When the flow 
is highly super sonic and super Alfvenic the interaction 
can be described using the well known dynamical fric- 
tion formalism of Chandrasekhar (1943) and Tremaine & 
Weinberg (1984). However, when the flow speed is below 
that of the fast magnetosonic speed but significantly ex- 
ceeds the slow magnetosonic speed, as in our case, the 
energy exchange rate is found to be reduced. This is par- 
ticularly the case when the relative velocity between gas 
and protoplanet is parallel or nearly parallel to the mag- 
netic field. A similar reduction would be expected to apply 
to the accretion rate by the protoplanet. 

Using the conservation of energy applied to the fluid 
perturbations in a more general context, we find that the 
energy exchange rate with the orbit is determined by ad- 



vected and wave energy fluxes at the system boundaries, 
confirming the importance of the existence of propagating 
waves and the boundary conditions for the determination 
of the energy exchange rate with the orbit. 

In section [4] we discuss our results, showing that, for 
parameters expected for protoplanetary disks, these im- 
ply negligible orbital evolution for protoplanets interior to 
magnetospheric cavities. This in turn implies that if entry 
into a magnetospheric cavity is responsible for a lack of 
further migration of a single protoplanet, the orbital ele- 
ments could subsequently only be affected by stellar tides. 
Finally we summarise our conclusions. 

2. Basic Equations 

The basic equations are the equations of ideal MHD in a 
frame rotating with uniform angular velocity e351», where, 
adopting cylindrical polar coordinates (r, <p, z) with ori- 
gin at the central star, e3 denotes the unit vector in the 
z direction and SI, is the angular velocity of the central 
magnetospheric region, taken to coincide with that of the 
central star. 

The equation of motion can be written in the form 

9v 

— + v • Vv + 2fi*e 3 x v = f, (1) 
ot 

where the force per unit mass is given by 



VP ^ (V x B) x B 
f = V$ + ^ —L , 

p VKp 

the induction equation in the form 
dB 



dt 



= V x (v x B) 



(2) 



(3) 



(4) 



and the continuity equation in the form 
|+V.(pv)=0. 

Here p is the density, P is the pressure, v is the velocity 
and the magnetic field is B. 

When the only source of gravity is the central star treated 
as a point mass, the combined gravitational and centrifu- 
gal potential is 



$ = — 



y/r 2 + 



with the central stellar mass being M» and G the gravita- 
tional constant. For the model considered here we adopt 
an isothermal equation of state so that P = pc 2 with c 
being the constant sound speed. We take this to be char- 
acteristic of the outer disk so that it will be much less than 
the orbital speed or a characteristic Alfven speed. 

3. Perturbation Due to an Orbiting Protoplanet 

3.1. Angular Momentum Transport and Equilibrium 
Stellar Spin 

We consider the perturbation of a magnetosphere interior 
to an accretion disk through which matter accretes at a 



J.C.B. Papaloizou: Protoplanet Magnetosphere Interactions 



3 



constant or slowly varying rate by an orbiting protoplanet. 
When a quasi-steady state has been reached, the disk is 
expected to be disrupted at a radius, Rd slightly interior to 
the corotation radius, ro where the stellar angular velocity 
coincides with the Keplerian angular velocity expected in 
the disk. This is given by 

r = (GM./SII) 1 ' 3 . (6) 

In the inner region the flow is expected to become mag- 
netically dominated and sub Alfv'cnic (Ghosh & Lamb 
1978). This general behaviour with the magnetospheric 
radius being ~ 0.85ro has been found in recent simula- 
tions (see Bouvier et al. 2006). It has also been estimated 
that the time required to attain such a quasi steady state 
is short compared to the lifetime of the disk (e.g. Cameron 
& Campbell 1993; Armitage & Clarke 1996). For a steady 
state flow and ^ imply that v and B are parallel. 

3.2. Local Analysis 

Finding the response of a general steady state magneto- 
spheric accretion flow to an orbiting protoplanet is a very 
difficult problem. Accordingly we consider possible sim- 
plifications. In this section we consider the calculation of 
the response using a local approximation. The idea behind 
this is that at any stage, the protoplanet moves relative 
to field lines along which the accretion flow moves, such 
that the speed of the protoplanet relative to the gas, de- 
noted by |U — v| below, may become comparable to the 
orbital speed (see Figure 1). The characteristic scale of 
interaction at which a protoplanet of mass, M p , would 
dominate the gas flow would then be set by the classical 
Bondi-Hoyle radius (e.g. Hoyle & Lyttleton 1939; Ruffert 
1996) r BH = 2GM p /\XJ - v| 2 . When the relative speed 
is characteristically orbital, tbh is smaller by a factor 
~ Mp/M* <C 1 than the radius at which the protoplanet 
orbits indicating a local interaction on a characteristic 
time scale that is smaller than the orbital time by a simi- 
lar factor. The smallness of teh/th ~ (M p /M*) 2 / 3 , with 
Th being the Hill radius also indicates that the interaction 
of the protoplanet with gas that is not bound to it is lin- 
ear even in the giant planet regime with M p /M* ~ 10~ 3 . 
The margins involved are such that the interaction is ex- 
pected to be local and linear even when the relative speed 
is significantly less than orbital. 

This is unlike the situation in protoplanetary disks 
where the small relative velocity between gas and proto- 
planet can make the interaction nonlinear at scales larger 
than the Hill radius resulting in gap formation (see e.g. 
Lin & Papaloizou 1993; Papaloizou & Terquem 2006; and 
references therein for more discussion). 

On the basis of the above discussion, we here adopt 
the assumption that the protoplanet magnetosphere inter- 
action is localised in the vicinity of the protoplanet and 
that the associated characteristic time scale is significantly 
shorter than the orbital period. Accordingly we adopt a 
local frame that can be taken to be a non rotating local 
Cartesian coordinate system with origin instantaneously 



at the centre of mass of the protoplanet and moving with 
a constant velocity equal to that of the rotational velocity 
the star would have if it extended to the radial location of 
the protoplanet (if one uses a local frame fixed in the orig- 
inal rotating frame and then neglects f2* in the equations 
of motion the same results are obtained). The response to 
the protoplanet is readily calculable and leads to a confir- 
mation of the conclusion that the changes induced in the 
orbit are associated with wave energy fluxes propagating 
away from the protoplanet. Furthermore when such waves 
are absent, such as when the relative speed is below all 
possible wave speeds, there is no energy change induced 
on the orbit. 

3.3. Linear Analysis 

The effect of a perturbing protoplanet of mass M p 
at position vector r = r or b is to produce an 
Eulerian perturbation to the gravitational potential $' = 
-GM p /,J\v-Y orb \ 2 + b 2 , where b is a softening param- 
eter which can be specified to ensure $' remains finite. 
An appropriate value is expected to be of the order of the 
radius of the planet. 

In the original rotating frame, the Lagrangian displace- 
ment £ induced by by this perturbing potential satisfies 
the equation (see eg. Frieman & Rotenberg 1960; Lynden- 
Bell & Ostriker 1967) 

D 2 £ D£ 

where the convective derivative operator 

< 8 » 

and the Lagrangian variation of the force per unit mass is 
F(€) = f + £ ■ Vf , with 

f = _|_ (yp - (V x B) x B 



+ 



p p* \ Ak 
(V x B') x B (V x B) x B' 



A'Kp 



A'Kp 



V$'. (9) 



The Eulerian perturbations to the magnetic field and den- 
sity are given by 



B' = V x (£ x B), 
and 



(10) 



(11) 



respectively. For an isothermal equation of state the 
Eulerian pressure perturbation P' = p'c 2 . Equations ([7] 
- 1111) thus provide a complete system for determining the 
Lagrangian displacement 

We make the usual local approximation that the steady 
state variables p, $ and P as well as v and B can be 
taken to be constant. In the local frame ( see Figure 1) 
the velocity of the protoplanet is taken to be U and we 
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Protoplanet 



Fig. 1. This illustration shows the protoplanet viewed in a frame rotating with the angular velocity of the star f2* here 
= f2, the latter being located at the origin of a Cartesian coordinate system (ei,e2,e3) defining the direction of the 
axes. The protoplanet moves in a circular orbit in the equatorial plane with velocity U. The field line passing through 
its location points out of the plane ultimately reaching the star. In this frame the gas moves along field lines with 
velocity v and the line pointing in the direction of the relative velocity v U makes an angle ip with that pointing 
in the direction of B, or equivalently Va, which also gives the direction of the z axis for the local frame in which the 
calculation is done. In that frame v — U lies in the (x, z) plane. 



may use equations ([7] - fTTT) with fi* 
then gives 



0.. Equation (0 



VP' 
P 



V$' 



(V x B') x B 

4-7T/9 



(12) 



which together with equations (fT0|) . (fTTj) and P' — —p'c 2 
enable the Lagrangian displacement induced by the per- 
turbing potential $' to be calculated. 

3.3.1. Fourier Transform of the Potential 

For this local model, we set the softening parameter 6 = 0, 
so that in the local coordinates, the perturbing potential 
may be written $' = — GM p /\r — Ut|, with r denoting the 
position vector. To solve (fl~2|) we take Fourier transforms of 
the perturbations. Thus for a general perturbation quan- 
tity Q, 



exp(ik • r)Q k (k)d 3 k, 

with the Fourier transform 
1 



(27T) 



exp(— ik • r)Q(r)d r. 



(14) 



The wavenumber k = (fc x , k yi k z ) and the integrals are 
taken over all coordinate and wavenumber space. For the 
perturbing potential one readily finds that 



GM P 

~27T 2 |k| 2 



exp(— ik • Ut). 



(15) 



3.3.2. Calculation of the Fourier Transform of the 
Displacement 

After taking the Fourier transforms of equation (|12p and 
equation (jTOJ) and noting that P k = p' w c 2 = — pc 2 £ k • k, 
one obtains a set of linear algebraic equations for the de- 
termination of the components of £ k in terms of $ k . The 
scalar quantity £ k k, is readily found after straightforward 
algebra to be given by 



^k 2 (. 2 -(k.V A ) 2 ) 



(yi - ( c 2 + V%)v 2 k 2 + (k • V A ) 2 c 2 k 2 ) 



(16) 



where v = (k ■ (v — U)) and Va = B/y^iTrp denotes a 
vector, with magnitude equal to the Alfven speed, and 
the direction of B. We set Va — |V A |. 

The gravitational force acting on the protoplanet as a 
result of the perturbation of the magnetosphere can then 
be found using 



(13) F orb = J P Wd 3 r. 



(17) 



We evaluate the rate of change of orbital energy that 
would be viewed in the fluid rest frame. In that frame the 
drag force would be zero if v = U and consequently the 
rate of change of orbital energy would be zero. Otherwise 
it is given by 

dE m 



dt 

dE or i, 



(U - v) ■ F orb or 



(18) 



dt 



= -(27r) 3 Pe / p(£ k ■ k)(k • (U - v))$-d 3 k. (19) 
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Here Re indicates that the real part of the integral is to 
be taken. 

In evaluating equation (|19|) we note that if there were 
no singularities in the expression for £ k • k given by equa- 
tion (p~6|) we would have no energy exchange or dE or i, /dt = 
0. Thus in order to obtain dE or b/dt ^ 0, there should be 
singularities in £ k ■ k. In fact singularities in the form of 
simple poles occur when the denominator on the right 
hand side of equation (|16[) vanishes or when 



(c 2 + v%yk 2 + (k • v A ) 2 c 2 k 2 



0. 



(20) 



This quartic gives four roots for v given by v — ±7"i 2 |k| , 
where r\ = (c 2 + V|)/2 + J(c 2 + Uf) 2 /4 - (k • V A ) 2 c 2 



andr 2 = (c 2 + Uf)/2- ^/(c 2 + V|) 2 /4 - (k • V A ) 2 c 2 , with 
k = k/|k|. 

The roots corresponding to r\ and r 2 , which without 
loss of generality are taken to be positive, correspond to 
fast and slow magnetosonic waves respectively which prop- 
agate with angular frequency v = k • (v — U) . When any 
such waves can be present, there is the possibility of an in- 
duced change of orbital energy. This is entirely consistent 
with the fact, noted in section [3.4.11 below that the rate of 
change of orbital energy has to be associated with a con- 
served energy flux propagating away from the protoplanet. 
In the local problem considered here, this flux can be car- 
ried in either fast or slow magnetosonic waves. We also 
comment that in order to produce orbital energy changes, 
as is apparent from equation (|17p. the waves have to be 
compressive or associated with density perturbations, thus 
incompressible Alfven waves are ineffective.. 

The singularities in £ k • k given by (fTB) can be dealt 
with by writing the denominator as a product of factors of 
the form (£±rj), for (i = 1, 2) v = f /|k| and then handling 
the inverse of each these by using the well known Landau 
prescription, thus 



1 



P 



1 



(21) 



{p±Ti) v(z/±ri). 

Here P denotes that the principal value is to be taken on 
integration and 5 is Dirac's <5 function. This prescription 
results from adding an infinitesimally small negative imag- 
inary part to v, which can be regarded as ensuring that 
the perturbing potential vanishes at t = — 00 thus impos- 
ing causality. Using the Landau prescription and noting 
that roots of opposite sign give equivalent contributions, 
we can find the rate of change of orbital energy using 
equation (|19[) in the form 

dEori 



dt 
where 



l^l^lkl 



r\r 2 (r{ - rf ) 



A(k)d 3 k, 



-n(r 2 2 -(k.V A ) 2 )^-r 2 ) 



(22) 



(23) 



To perform the above integral we first remark that because 
r\ , r 2 and v are functions only of the unit wavenumber vec- 
tor, if polar coordinates in wavenumber space are adopted, 



the integrand easily factors into radial and angular parts, 
thus 



dE, 



orb 



dt 



-2(GM p fp ^ — 



z>A(k) 



r\r 2 (r{ - r 2 ) 



2 dCl, (24) 



where dfl denotes integration over the solid angle in 
wavenumber space and k max , k min are the magnitudes of 
upper and lower wavenumber cut offs that necessarily have 
to be used in local calculations of gravitational interac- 
tions because of the long and short range form of the 1 jr 
potential (see eg. Tremaine & Weinberg 1984 in the con- 
text of stellar dynamics and Papaloizou 2002 in the con- 
text of thin hydrodynamic disks). Physically k^ ax should 
be taken as the smallest effective length scale in the prob- 
lem and the largest. For a giant protoplanet in an 
inner magnetosphere, the ratio of the radius of the planet 
to a characteristic scale of the magnetosphere of 0.1 AU, 
r-j 10~ 2 , accordingly one expects k max /k m i n ~ 100. 

To perform the angular integral in eauation (|24p . we 
specify the local frame such that the z axis points 
along the direction of B and thus V A and v. The 
direction of v — U is taken to be in the (x, z) plane 
making an angle ip with the z axis (see Figure 1). 
Using spherical polar angles 9, </> to perform the in- 
tegral over the solid angle, we note that because 
k • V A /Vyt = cos#, ri,r 2 are functions only of 9. Given 
that v = |v — U| (cos 9 cos ip + sin 9 sin ip cos <fi ) , the 
delta functions readily enable the integration over <f> and 
so we obtain an expression requiring only an integration 
over \x = cos 9 in the form 

dE orb Att(GM p ) 2 P 



dt 



Iv-U| 



In 



V, 



where T> = (I\,2 + ^2,1)/^, with 

(r 2 - Vif)dn 



I 



1,2 



(rf - r 2 )^J (1 - f 2 ) sin 2 ^ - (ji - fi cos i>f 



(25) 



(26) 



where (/1, f 2 ) = (r 1; r 2 )/|v — U|. In the above integral the 
domain of integration is either (—1,1) or (s 1 ,s 2 ) ; where 
si and S2 are values of fi for which the square root in (|26|) 
is zero should they exist. 

In general the dependence of /1 and f 2 on p through 
r\ and r 2 makes the above integral intractable. However, 
it can be found in the limit of low c/Va that we consider 
here. Then r± =Va, and r 2 = cfi. 

We first note that In the super Alfvenic case where 
both /1 and f 2 are very small, one obtains T> = 1. This 
case of course applies to the limit of zero wave propaga- 
tion speeds which corresponds to the dynamical friction 
calculation of Chandrasekhar (1943) and it can also be 
derived from the calculation of the linear response of col- 
lisionlcss particles by Tremaine & Weinberg (1984). We 
obtain D < 1 so that then this quantity measures the 
factor by which the dynamical friction is reduced when 
compared to the estimate of Chandrasekhar (1943). 

To emphasise that this can be a large reduction, we 
note that in the sub Alfvenic case when both /1 and f 2 
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exceed unity, there is no wavelike response or domain of 
integration for (-Zi^, ^2,1) so that T> = 0. 

When fx is large but ji very small so that slow waves 
exist but fast waves do not, we have I12 = and l2,i = 
(tt sin t/>)/2 so that D = (sin 2 ^)/2. 

Thus in this case the rate of change of orbital energy de- 
pends on the orientation of the magnetic field, being zero 
when t/> = 0. In this limit U, v and B are all parallel. 

There is no dynamical friction in that case because de- 
flection of a fluid clement towards the protoplanet requires 
motion perpendicular to the magnetic field lines. 
In contrast when ip = tt/2 deflection towards the proto- 
planet can occur with the motion remaining parallel to 
the field lines. Then the interaction is uninhibited. 
Finally when /1 < 1 and / 2 = 0, so that both waves exist 
but the ratio of propagation speeds is large, we obtain 
V = 1 — / 2 (cos 2 ip — (sin 2 ip)/2). In this case the most 



effective reduction of the dynamical friction occurs when 
ip = 0, then V = 1- ff. 

3.4. General Scalings for the Response and Torque 

3.4.1. Conservation of Wave Action and Angular 
Momentum Evolution 



source Sa is localised in space, the rate of change of orbital 
energy is associated with a conserved energy flux propa- 
gating away to large distances from the protoplanet. The 
total effect of the protoplanet is obtained by summing the 
independent contributions from different m. 

When the unperturbed configuration is axisymmetric 
and the perturbing potential $' depends on tp and t in the 
combination uj p t — ip, the energy fluxes convert to fluxes 
of the angular momentum component along the symmetry 
axis by dividing by the pattern speed. 

Then the total torque acting on the fluid can be writ- 
ten as T = dJfiuid/dt = Wp dE fi U id/ dt. When, as for 
the situations considered here, the wave fluxes produce an 
energy and angular momentum loss from the system, the 
orbit decays or dE or i,/dt = —dEfi U id/dt < 0. 

We comment that the fact that the energy changes in 
the orbit can be measured through energy fluxes at dis- 
tant boundaries generalises a corresponding result for hy- 
drodynamic disk forcing (see Papaloizou & Terquem 2006) 
to the general MHD case. In particular if there are no ex- 
cited waves or advected disturbances, there are no induced 
changes to the orbit. 

3.4.2. Torque Scaling 



We point out that in the inviscid case considered here, it 
is possible to derive a conservation law from equations ([7J 
- [TT1) in the form 



^ + V • ( P av + F A ) = S A , 



(27) 



where in Cartesian coordinates with r = (x%, x%, X3), v = 
(v\, V2,V3), Im denoting that the imaginary part is to be 
taken and using the summation convention 

dt 



pA = Im\ pi 



F A = Im 



UL 

B'*x(£ x B) 



(e 3 x£)-£ 



d(pvjt) 



dxj 



4-7T 



V x B ■ (C x g)B 

8^ 



and 



The scaling of the rate of energy change given by equation 
([251) . that was obtained from the local analysis, with the 
physical parameters of the problem can be obtained quite 
generally. The length scale appropriate to the magneto- 
sphere and the orbit is expeced to be ro . The unit of time 
is f2~ , and the pattern speed is expected to be u p ~ f2*, 
giving a characteristic relative velocity rofi*. This veloc- 
ity would be expected to be characteristic of both the flow 
velocity along field lines and the relative velocity between 
(28) ^ ne or biting protoplanet and the magnetosphere. For a 
characteristic density in the neighbourhood of the proto- 
planet, p, we find natural scalings from equations ([71- [TTj) 
(29)and (j3Tj) for the response displacement and torque such 
that | oc r M p /M* and T oc p(GM p ) 2 / (Q*r ). Thus we 
write 



-Im (pC ■ V$') 



(30) 



T 



This conservation law is written for a general complex 
forcing potential Here we apply it to the situation where 
the forcing can be written as the real part of a potential 
Re ($') , where $' = g m (r) ex^>{im{uj p t — tp)), g m is a com- 
plex function of r, m the azimuthal mode number and tu p 
is the pattern speed, being the difference between the an- 
gular velocity of the orbiting protoplanet and f2». In this 
case, provided the forcing potential vanishes at the fluid 
boundaries, the rate of increase of the energy of the fluid 
is 

Re { p>)Re(^)d*r = ^[s A d*r.(31) 



-1 dE or 

>n dt 



fp(GM p ) 2 /(nlr Q ), 



(32) 



where / includes dependence on softening or the small 
scale cut off as well as the existence of propagat- 
ing waves. This is equivalent to quation (|25[) if / = 
4irr nl In (fc 

max /kmin) (|v - U\uJ p ) 



dEfi U id 



dt 



dt 



4. Discussion 

4.1. Orbital Evolution Timescale 

From equation (f3"2"]> a characteristic rate of evolution of 
protoplanet in circular orbit at radius r or i, can be esti- 
mated from 



The quantities muj p pA/'2,muj p Fj^/2 may be regarded as 
the energy density and the non advected energy flux as- 
sociated with the forcing. In particular when the forcing 



torb ~ 



dE or i> j dt 



E, 



orb 



2pG 2 M p r orb 




(33) 
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For the magnetically dominated region (see eg. Bouvier 
et al. 2006) we set p = M/(47irgfi»). Then for M» = 
1Mq,M p = 10~ 3 M*, r or j = 0.5r giving a location 
slightly interior to the 2:1 commensurability with fi*, 
and a characteristic accretion rate for a protostellar disk 
M = 3 x lO" 8 M y^ 1 ( e.g. Muzerolle et al. 2003 ), we ob- 
tain t rb ~2x 10 11 y. Thus the expected evolution of such 
a protoplanet orbit is expected to be very small over a 
characteristic protostellar disk lifetime ~ I0 7 y. Note too 
that this result is not changed if the orbital decay rate 
is enhanced by a factor / ~ 50, corresponding to use of 
equation ([25]) with T> = 1 corresponding to effective wave 
propagation, |v - U| = (r £l*)/2 and k max /k mm ~ 100. 
It also indicates that no additional mechanism, such as 
special torques acting near the inner disk edge (Masset et 
al. 2006) is needed to halt the inward migration of proto 
giant planets. 

4.1.1. Protoplanet Accretion 

It may also be argued that the expected accretion onto 
the protoplanet is negligible while it is inside the magne- 
tosphcre. In order that accretion can take place an amount 
of energy comparable to the orbital binding energy must 
be dissipated in order for material to become bound to the 
protoplanet. Thus an estimate of the accretion rate onto 
the protoplanet is given by 
1 dE oro 



M p ~ M, 



E, 



orb 



dt 



(34) 



which implies that the accretion time scale is the same as 
that for orbital evolution. 



4.1.2. Other Effects 

We here consider other effects that may lead to orbital 
decay of the protoplanet. First we derive an effective 
drag coefficient characterising the dynamical friction act- 
ing through the gravitational torques calculated above. 
To do this we write 

,1E - U V^/^IJ v:i; v;-i; 



dt 2 ~...-^~ - (35) 

( Landau & Lifshitz 1993), where the expression is ap- 
plied locally and Cp, is the drag coefficient, f Adopting 
tq/R p = 200, we find that the larger estimate of t or b given 
above corresponds to Cd ~ 0.1 while the smaller estimate 
corresponds to Cd ~ 5. These results imply that the or- 
bital decay rate is, to within an order of magnitude or so, 
comparable to the non gravitational effects arising from 
the protoplanet acting as an obstacle to the flow. This is 
in contrast to the situation that occurs for protoplanets 
in circular orbits in a thin disk (Lin & Papaloizou 1993). 

4.1.3. Magnetic Coupling with the Central Star 

Orbital energy decay can also result from the direct inter- 
action of a conducting protoplanet with the stellar mag- 
netic field when this varies around the orbit. When the 



stellar magnetic field is non axisymmetric and the proto- 
planet has a non zero resistivity , there will be periodic 
flux penetration and dissipation leading to an orbital en- 
ergy loss. In the discussion below we neglect any internal 
protoplanetary magnetic field. 

The rate of change of energy can be estimated to be 
(see e.g. Joss, Katz & Rappaport 1979 for a discussion of 
magnetic binary systems) 



dE n 



dt 



-R 2 p B^AR p \ 



(36) 



where ARp = ^Jr)]{w v — v^/r) is a measure of the depth 
of external flux penetration per relative orbit, with v v 
being the azimuthal component of the gas velocity, B\ \ is 
a typical component of B parallel to the surface of the 
protoplanet, which we shall take to be |B| and r\ is its 
magnetic diffusivity. We may also write 



dE n 



3?r ^3 D 2 Ai V 



dt 



-pVlK- 



vjr) 



V A 



which leads to an effective drag coefficient given by 



Co = 



SVjrj^p - v v /r) 
|U- v|(U- v) -U| 



R v \ ( AR Z 



(37) 



(38) 



Although the flow is sub Alfvenic, the above is of or- 
der unity apart from the two factors in brackets. Of these 
Rp/r ~ 10~ 2 and the last factor which represents the ra- 
tio of the skin depth to radius is also small. For example 
from Zhang, Jones & Chen (1996) the magnetic diffusiv- 
ity of Jupiter can be estimated to be rj = 3 x 10 3 cm 2 s _1 . 
For a relative orbital period of 10 days, this leads to 
ARp/Rp ~ 10~ 5 and the effective C D ~ 10~ 6 . Thus the 
effect of this type of magnetic interaction with the central 
star is not likely to be significant. 

4.2. Conclusions 

In this paper we performed a local calculation of the proto- 
planet magnetosphere interaction derived. The existence 
of propagating waves was necessary for there to be a non 
zero energy exchange rate. This is also expected from very 
general considerations of energy and wave action con- 
servation. For high relative speeds between the gas and 
protoplanet, the interaction was found to be identical to 
that obtained from the dynamical friction formalism of 
Chandrasekhar (1943). But when, as expected here, the 
relative flow speed is below the fast magnetosonic speed 
but exceeds the slow magnetosonic speed, the interaction 
strength was found to be reduced, especially when the rel- 
ative velocity between gas and protoplanet was parallel or 
nearly parallel to the magnetic field. The inhibition of the 
ability of the protoplanet to disturb the flow in such cases 
would also be expected to apply to the gas accretion rate 
onto it. 

For parameters expected for protoplanetary disks our 
calculation indicates negligible orbital evolution for proto- 
planets interior to magnetospheric cavities through which 
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accretion is taking place during expected protoplanetary 
disk lifetimes. This is the situation when the protoplan- 
ets are far enough away from the inner disk edge so that 
interaction with the disk is negligible and accordingly it 
is not necessary to invoke special torques associated with 
the disk edge to halt or reverse migration. In this regard 
we comment that although there have been, of necessity, 
many simplifying assumptions made in order to carry out 
the analysis presented here, the expected protoplanet or- 
bital evolution in the magnetosphere fails to be significant 
by a wide margin. 
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